February 2, 2008 5:34 WSPC/Lecture Notes Series: 9in x 6in larson-chapter 



o 
o 

< 

OO 

^S| ■ Unitary Systems, Wavelet Sets, and Operator-Theoretic 

Interpolation of Wavelets and Frames 

< 

-! ' David R. Larson 

>-Q ' 
■ 4—* ' Department of Mathematics 

' Texas A&M University 

College Station, TX 77845 

E-mail: larson@math.tamu.edu 

. 'iL ' A wavelet is a special case of a vector in a separable Hilbert space that 

I I generates a basis under the action of a collection, or system, of uni- 

VQ ■ tary operators. We will describe the operator-interpolation approach to 

^^ ' wavelet theory using the local commutant of a system. This is really an 

("^ , abstract application of the theory of operator algebras to wavelet theory. 

V^ . The concrete applications of this method include results obtained using 

("V^ ' specially constructed families of wavelet sets. A frame is a sequence of 

fH , vectors in a Hilbert space which is a compression of a basis for a larger 

space. This is not the usual definition in the frame literature, but it 

j is easily equivalent to the usual definition. Because of this compression 

P^ . relationship between frames and bases, the unitary system approach to 

wavelets (and more generally: wandering vectors) is perfectly adaptable 

to frame theory. The use of the local commutant is along the same lines 

^^ ■ as in the wavelet theory. Finally, we discuss constructions of frames with 

$_( ' special properties using targeted decompositions of positive operators, 

5^ , and related problems. 

1. Introduction 

This is a write-up of of a tutorial series of three talks which I gave as part 
of the "Workshop on Functional and Harmonic Analyses of Wavelets and 
Frames" held August 4-7, 2004 at the National University of Singapore. I 
will first give the titles and abstracts essentially as they appeared in the 
workshop schedule. I will say that the actual style of write-up of these notes 
will be structured a bit differently, but only in that more than three sections 
will be given, and subsections indicated, to (hopefully) improve expositional 
quality. 



B 
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1.1. Talks and Abstracts 

(a) " Unitary Systems and Wavelet Sets" : A wavelet is a special case of a 
vector in a separable Hilbert space that generates a basis under the action 
of a collection, or "system" , of unitary operators defined in terms of trans- 
lation and dilation operations. This approach to wavelet theory goes back, 
in particular, to earlier work of Goodman, Lee and Tang [25] in the con- 
text of multiresolution analysis. We will begin by describing the operator- 
interpolation approach to wavelet theory using the local commutant of a 
system that was worked out by the speaker and his collaborators a few 
years ago. This is really an abstract application of the theory of operator 
algebras, mainly von Neumann algebras, to wavelet theory. The concrete 
applications of operator-interpolation to wavelet theory include results ob- 
tained using specially constructed families of wavelet sets. In fact X. Dai 
and the speaker had originally developed our theory of wavelet sets [11] 
specifically to take advantage of their natural and elegant relationships 
with these wavelet unitary systems. We will also discuss some new results 
and open questions. 

(b) " Unitary Systems and Frames" : A frame is a sequence of vectors in a 
Hilbert space which is a compression of a basis for a larger space. (This is not 
the usual definition in the frame literature, but it is equivalent to the usual 
definition. In this spirit, the usual "inequality" definition can be thought 
of as an abstract characterization of a compression of a basis.) Because of 
this compression relationship between frames and bases, the unitary system 
approach to wavelets (and more generally: wandering vectors) is perfectly 
adaptable to frame theory. This idea was developed into a theory a few 
years ago by D. Han and the speaker [33]. The use of the local commutant 
is along the same lines. 

(c) "Decompositions of Operators and Operator- Valued Frames": We 
will discuss some joint work with K. Kornclson and others on construction 
of frames with targeted properties [16, 42]. These are related to targeted 
decompositions of positive operators. 



1.2. Some Background 

It might be appropriate to give some comments of a personal-historical 
nature, before continuing with the technical aspects. My particular point 
of view on "wavelet theory", which was developed jointly with my good 
friend and colleague Xingde Dai, began in the summer of 1992. Before then, 
I was strictly an operator theorist. I had heard my approximation theory 
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colleagues and friends at Texas A&M University talk about wavelets and 
frames, and Xingde had frequently mentioned these topics to me when he 
was finishing his Ph.D. at A&M (he was my student, graduating in 1990, 
with a thesis [9] on the subject of nest algebras). But it was this meeting 
of minds we had in June of 1992 that was the turning point for me. We 
formulated an approach to wavelet theory (and as it turned out ultimately, 
to frame- wavelet theory as well) that I felt we could " really understand" as 
operator algebraists. This was the abstract unitary system approach. Dai 
knew the unitary operator approach to multiresolution analysis that had 
been recently (at that time) published by Goodman, Lee and Tang 25], and 
he suggested to me that we should try to go further with these ideas in an 
attempt to get a some type of tractable classification of all wavelets. We 
went, in fact, in some completely different directions. The first paper that 
came out of this was the AMS Memoir [11] with Dai. The second paper was 
our paper [12] with Dai and Speegle, which proved the existence of single 
wavelets in higher dimensions, for arbitrary expansive dilations. After that, 
several papers followed including the AMS Memoir [33] with Deguang Han, 
and the Wutam Consortium paper [52], as well as the papers [10, 13, 26, 
27, 28, 30, 34, 43, 44] by my students Dai, Gu, Han and Lu, and their 
collaborators, and the papers [2, 39, 40] with my colleagues Azoff, lonascu, 
and Pearcy. 

The paper [11] with Dai mentioned above, which was published in 1998, 
culminated about two years of work on this topic by the authors. It con- 
tained our entire operator-theoretic approach to wavelet theory, and was 
completed in December 1994. This work, while theoretical, had much hands- 
on experimentation in its development, and resulted in certain theorems we 
were able to prove concerning constructions of new families of wavelets. In 
order to conduct successful experiments with our operator techniques, we 
needed a supply of easily computable test wavelets: that is, wavelets which 
were very amenable to paper and pencil computations. We discovered that 
certain sets, we called wavelet sets, existed in abundance, and we computed 
many concrete examples of them along the way toward proving our results 
of [11]. Several of these were given as examples in [11, Example 4.5, items 
(i) -^ (xi)]. Some of these are given in Section 2.6.1 of the present article. 

Most of our work in [11] on the local commutant and the theory of 
wavelet sets was accomplished in the two-month period July-September 
1992. The first time Dai and I used the terms wavelet set and local com- 
mutant, as well as the first time we discussed what we referred to as the 
connectedness problem for wavelets, was in a talk in a Special Session on 
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Operator Algebras in the October 1992 AMS Sectional Meeting in Dayton. 
Along the way a graduate student at A&M, Darrin Speegle, who was 
enrolled in a seminar course of Larson on the manuscript of [11], answered 
an open question Larson gave out in class by proving that the set of all 
wavelet sets for a given wavelet system is connected in the symmetric differ- 
ence metric on the class of measurable sets of finite measure. That resulted 
in a paper [50] which became part of his thesis (which was directed by 
William Johnson of Texas A&M), and Speegle subsequently joined forces 
with Dai and Larson [12] to prove that wavelet sets (and indeed, wavelets) 
exist in much greater generality than the prevailing folklore dictated. We 
received some attention for our work, and especially we thank Guido Weiss 
and John Benedetto for recognizing our work. This led to a flurry of papers 
by a number of authors, notably [4, 6, 13], and also led to the paper [52] 
by the Wutam Consortium, which was a group led by Guido Weiss and 
Larson, consisting of 14 researchers-students and postdocs of Weiss and 
Larson-based at Washington University and Texas A&M University, for 
the purpose of doing basic research on wavelet theory. 

1.2.1. Interpolation 

The main point of the operator-theoretic interpolation of wavelets (and 
frames) that Dai and I developed is that new wavelets can be obtained as 
linear combinations of known ones using coefficients which are not neces- 
sarily scalars but can be taken to be operators (in fact, Fourier multipliers) 
in a certain class. The ideas involved in this, and the essential computa- 
tions, all extend naturally to more general unitary systems and wandering 
vectors^ and I think that much of the theory is best-put in this abstract set- 
ting because clarity is enhanced, and because many of the methods work 
for more involved systems that are important to applied harmonic analysis, 
such as Gabor and generalized Gabor systems, and various types of frame 
unitary systems. 

1.2.2. Some Basic Terminology 

This article will concern bounded linear operators on separable Hilbert 
spaces. The set of all bounded linear operators on a Hilbert space H will be 
denoted by B{H). By a bilateral shift U on H we mean a unitary operator U 
for which there exists a closed linear subspace E G H with the property that 
the family of subspaces {U"E : n G Z} are orthogonal and give a direct- 
sum decomposition of H. The subspace E is called a complete wandering 
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subspace for U. The multiplicity of U is defined to be the dimension of E. 

The strong operator topology on B{H is the topology of pointwise con- 
vergence, and the weak operator topology is the weakest topology such that 
the vector functionals cux.y on B{H) defined by A k^ {Ax,y), A e B{H), 
x,y ^ H, are all continuous. An algebra of operators is a linear subspace 
of B{H) which is closed under multiplication. An operator algebra is an 
algebra of operators which is norm-closed. A subset S C B{H) is called 
selfadjoint if whenever A (z S then also A* € S. A C*- algebra is a self- 
adjoint operator algebra. A von Neumann algebra is a C*-algebra which is 
closed in the weak operator topology. For a unital operator algebra, it is 
well known that being closed in the weak operator topology is equivalent 
to being closed in the closed in the strong operator topology. 

The commutant of a set S of operators in B{H) is the family of all 
operators in B{H) that commute with every operator in S. It is closed 
under addition and multiplication, so is an algebra. And it is clearly closed 
in both the weak operator topology and the strong operator topology. We 
use the standard prime notation for the commutant. So the commutant of 
a subset S C B{H) is denoted: S' := {A e B{H) : AS = SA, S e S}. 

The commutant of a selfadjoint set of operators is clearly a von Neu- 
mann algebra. Moreover, by a famous theorem of Fuglede every operator 
which commutes with a normal operator N also commutes with its adjoint 
N* , and hence the commutant of any set of normal operators is also a von 
Neumann algebra. So, of particular relevance to this work, the commutant 
of any set of unitary operators is a von Neumann algebra. 

One of the main tools in this work is the local commutant of a system 
of unitary operators. (See section 2.4.) This is a natural generalization of 
the commutant of the system, and like the commutant it is a linear space 
of operators which is closed in the weak and the strong operator topologies, 
but unlike the commutant it is usually not selfadjoint, and is usually not 
closed under multiplication. It contains the commutant of the system, but 
can be much larger than the commutant. The local commutant of a wavelet 
unitary system captures all the information about the wavelet system in an 
essential way, and this gives the flavor of our approach to the subject. 

If [/ is a unitary operator and A is an operator algebra, then U is said 
to normalize A if U* ■ A-U = A . In the most interesting cases of operator- 
theoretic interpolation: that is, those cases that yield the strongest struc- 
tural results, the relevant unitaries in the local commutant of the system 
normalize the commutant of the system. 
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2. Unitary Systems and Wavelet Sets 

We define a unitary system to be simply a collection of unitary operators 
U acting on a Hilbert space H which contains the identity operator. The 
interesting unitary systems all have additional structural properties of var- 
ious types. We will say that a vector ^ G H is wandering for U if the set 



UtP := {U^j :U eU} (1) 

is an orthonormal set, and we will call ip a complete wandering vector for 
U if U')p spans H. This (abstract) point of view can be useful. Write 

W{U) 

for the set of complete wandering vectors for U. 

2.1. The One- Dimensional Wavelet System 

For simplicity of presentation, much of the work in this article will deal 
with one-dimensional wavelets, and in particular, the dyadic case. The other 
cases: non-dyadic and in higher dimensions, are well-described in the liter- 
ature and are at least notationally more complicated. 

2.1.1. Dyadic Wavelets 

A dyadic orthonormal wavelet in one dimension is a unit vector ip e 
L^(]R, /i), with /i Lebesgue measure, with the property that the set 

{2tV(2"t-0 :«,'eZ} (2) 

of all integral translates of ip followed by dilations by arbitrary integral 
powers of 2, is an orthonormal basis for L^(R, /i). The term dyadic refers 
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to the dilation factor "2". The term mother wavelet is also used in the 
literature for -tp. Then the functions 

are called elements of the wavelet basis generated by the " mother" . The 
functions ipn,i will not themselves be mother wavelets unless n = 0. 

Let T and D be the translation (by 1) and dilation (by 2) unitary 
operators in B{L^{R) given by (T/)(t) = f{t - 1) and {Df){t) = V2f{2t). 
Then 

25V(2"t-0 == (i:'"T'v)(t) 

for all n, I G Z. Operator-theoretically, the operators T, D are bilateral shifts 
of infinite multiplicity. It is obvious that L^{[0, 1]), considered as a subspace 
of L^(M), is a complete wandering subspace for T, and that L^([— 2, —1] U 
[1,2]) is a complete wandering subspace for D. 

2.1.2. The Dyadic Unitary System 

Let Ud,t be the unitary system defined by 

Ud,t ^ {D'^T^ -.nJeZ} (3) 

where D and T are the operators defined above. Then ip is a dyadic or- 
thonormal wavelet if and only if -0 is a complete wandering vector for the 
unitary system Ud,t- This was our original motivation for developing the 
abstract unitary system theory. Write 

WiD, T) := W{Ud^t) (4) 

to denote the set of all dyadic orthonormal wavelets in one dimension. 

An abstract interpretation is that, since D is a bilateral shift it has 
(many) complete wandering subspaces, and a wavelet for the system is a 
vector -0 whose translation space (that is, the closed linear span of {T^ : 
fc G Z} is a complete wandering subspace for D. Hence tp must generate 
an orthonormal basis for the entire Hilbert space under the action of the 
unitary system. 

2.1.3. Non-Dyadic Wavelets in One Dimension 

In one dimension, there are non-dyadic orthonormal wavelets: i.e. wavelets 
for all possible dilation factors besides 2 (the dyadic case). We said "possi- 
ble" , because the scales {0,1,-1} are excluded as scales because the dilation 
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operators they would introduce are not bilateral shifts. All other real num- 
bers for scales yield wavelet theories. In [11, Example 4.5 (x)] a family of 
examples is given of three-interval wavelet sets (and hence wavelets) for all 
scales d> 2, and it was noted there that such a family also exists for dila- 
tion factors 1 < d < 2. There is some recent (yet unpublished) work that 
has been done, by REU students and mentors, building on this, classifying 
finite-interval wavelet sets for all possible real (positive and negative scale 
factors). I mentioned this work, in passing, in my talk. 

2.2. N dimensions 

2.2.1. The Expansive-Dilation Case 

Let 1 < TO < oo, and let ^4 be an n x n real matrix which is expansive 
(equivalently, all (complex) eigenvalues have modulus > 1). By a dilation - 
A regular-translation orthonormal wavelet we mean a function ij) G -L^(M") 
such that 

{|dei(A)|^^(A"t - {hM. -, InY -.nJeZ} (5) 

where t = (ii, ..., t„)*, is an orthonormal basis for L^(R";to,). (Here m is 
product Lebesgue measure, and the superscript "t" means transpose.) 

If A € M„ (K) is invertible (so in particular if A is expansive) , then it is 
very easy to verify that the operator defined by 

(DAmt) = \detA\if{At) (6) 

for / e L^(R"), i e K", is unitary. For 1 < i < n, let T^ be the unitary 
operator determined by translation by 1 in the i*^ coordinate direction. 
The set (5) above is then 

{D'XTi'---T!^"^P:k,keZ} (7) 

If the dilation matrix A is expansive, but the translations are along 
some oblique lattice, then there is an invertible real n x n matrix T such 
that conjugation with Dt takes the entire wavelet system to a regular- 
translation expansive-dilation matrix. This is easily worked out, and was 
shown in detail in [39] in the context of working out a complete theory 
of unitary equivalence of wavelet systems. Hence the wavelet theories are 
equivalent. 

2.2.2. The Non-Expansive Dilation Case 

Much work has been accomplished concerning the existance of wavelets for 
dilation matrices A which are not expansive. Some of the original work was 
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accomplished in the Ph.D. theses of Q. Gu and D. Specglc, when they were 
together finishing up at Texas A&M. Some significant additional work was 
accomplished by Speegle in [49], and also by others. In [39], with lonascu 
and Pearcy we proved that if an nxn real invertible matrix A is not similar 
(in the nxn complex matrices) to a unitary matrix, then the correspond- 
ing dilation operator Da is in fact a bilateral shift of infinite multiplicity. 
If a dilation matrix were to admit any type of wavelet (or frame-wavelet) 
theory, then it is well-known that a necessary condition would be that the 
corresponding dilation operator would have to be a bilateral shift of infinite 
multiplicity. I am happy to report that in very recent work [45], with E. 
Schulz, D. Speegle, and K. Taylor, we have succeeded in showing that this 
minimal condition is in fact sufficient: such a matrix, with regular transla- 
tion lattice, admits a (perhaps infinite) tuple of functions, which collectively 
generates a frame-wavelet under the action of this unitary system. 

2.3. Abstract Systems 

2.3.1. Restrictions on Wandering Vectors 

We note that most unitary systems U do not have complete wandering 
vectors. For W{U) to be nonempty, the set U must be very special. It must 
be countable if it acts separably (i.e. on a separable Hilbert space), and it 
must be discrete in the strong operator topology because if U,V Cz U and 
if a; is a wandering vector for U then 

\\U-V\\ > \\Ux-Vx\\ =V2 

Certain other properties are forced on U by the presence of a wandering 
vector. One purpose of [11] was to study such properties. Indeed, it was a 
matter of some surprise to us to discover that such a theory is viable even 
in some considerable generality. For perspective, it is useful to note that 
while Ud,t has complete wandering vectors, the reversed system 

i^T,D = {T^D"^ :nJeZ} 

fails to have a complete wandering vector. (A proof of this was given in the 
introduction to [11].) 

2.3.2. Group Systems 

An example which is important to the theory is the following: let G be an 
arbitrary countable group, and let H — P{G). Let tt be the (left) regu- 
lar representation of G on H. Then every element of G gives a complete 
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wandering vector for the unitary system 

U:=n{G). 

{li h G G it is clear that the vector A/j G l'^{G), which is defined to have 1 in 
the h position and elsewhere, is in W{U).) If a unitary system is a group, 
and if it has a complete wandering vector, it is not hard to show that it is 
unitarily equivalent to this example. 

2.4. The Local Commutant 

2.4.1. The Local Commutant oj the System,UD,T 

Computational aspects of operator theory can be introduced into the 
wavelet framework in an elementary way. Here is the way we originally 
did it: Fix a wavelet tp and consider the set of all operators S E B{L'^{M.)) 
which comm,ute with the action of dilation and translation on ip. That is, 
require 

(5V')(2"t-0 = 5'(V'(2"t-/)) (8) 

or equivalently 

D^T^Sip = SD'^T^-iP (9) 

for all n, Z G Z. Call this the local commutant of the wavelet system IAd.t 
at the vector ip. (In our first preliminary writings and talks we called it 
the point commutant of the system . Formally, the local commutant of the 
dyadic wavelet system on L^(R) is: 

C-4,{Ud,t) :- {S e B{L^{R)) : (SD''T^ - D''T'S)i^ ^ 0,Vn,? e Z} (10) 

This is a linear subspace of B{H) which is closed in the strong operator 
topology, and in the weak operator topology, and it clearly contains the 
commutant of {D,T}. 

A motivating example is that if rj is any other wavelet, let V := V^ be 
the unitary (we call it the interpolation unitary) that takes the basis -ipn^i 
to the basis rjnj. That is, V^n,i — Vn.i for all n,l E Z,. Then 77 = Vip, so 
VD^-T^ijj = D'^T^Sij hence V e C^{Ud,t)- 

In the case of a pair of complete wandering vectors ip, rj for a general 
unitary system hi, we will use the same notation V2 for the unitary that 
takes the vector Uip to Urj for all U €U. 

This simple-minded idea is reversible, so for every unitary V in C^{Ud,t) 
the vector V^ is a wavelet. This correspondence between unitaries in 
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C^{D,T) and dyadic orthonormal wavelets is one-to-one and onto (see 
Proposition 1.) This turns out to be useful, because it leads to some new 
formulas relating to decomposition and factorization results for wavelets, 
making use of the linear and multiplicative properties oi C^{D,T). 

It turns out (a proof is required) that the entire local commutant of the 
system Ud,t at a wavelet ^ is not closed under multiplication, but it also 
turns out (also via a proof) that for most (and perhaps all) wavelets ip the 
local commutant at ip contains many noncommutative operator algebras (in 
fact von Neumann algebras) as subsets, and their unitary groups parame- 
terize norm-arcwise-connected families of wavelets. Moreover, C^{D,T) is 
closed under left multiplication by the commutant {D, T}' , which turns 
out to be an abelian nonatomic von Neumann algebra. The fact that 
C^{D,T) is a left module under {D,T}' leads to a method of obtaining 
new wavelets from old, and of obtaining connectedness results for wavelets, 
which we called operator-theoretic interpolation of wavelets in [DL], (or 
simply operator-interpolation). 

2.4.2. The Local Commutant of an Abstract Unitary System 

More generally, let S C B{H) be a set of operators, where iJ is a separable 
Hilbert space, and let x G iJ be a nonzero vector, and formally define the 
local commutant of 5 at x by 

C^{S) := {A e B{H) : {AS - SA)x = 0,SeS} 

As in the wavelet case, this is a weakly and strongly closed linear sub- 
space of B{H) which contains the commutant S' of S. If x is cyclic for S in 
the sense that span(5x) is dense in H, then x separates Cx{S) in the sense 
that for S E Cx{S), we have S'x = iff a; = 0. Indeed, if A e Cx{S) and if 
Ax = 0, then for any S £ S we have ASx = SAx = 0, so ASx =- 0, and 
hence A = 0. 

If A e Cx{S) and B e 5', let C = BA. Then for all S* G 5, 

(CS - SC)x ^ B{AS)x - {SB) Ax ^ B{SA)x - {BS)Ax = 

because ASx = SAx since A e Cx{S), and SB = BS since B G S' . Hence 
Cx{S) is closed under left multiplication by operators in S' . That is, Cx{S 
is a left module over S' . 

It is interesting that, if in addition 5 is a multiplicative semigroup, 
then in fact Cx{S) is identical with the commutant S' so in this case the 
commutant is not a new structure. To see this, suppose A G Cx{S). Then 
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for each S,T E S we have ST E S, and so 

AS(Tx) = {ST)Ax = S{ATx) ^ (S)Tx 
So since T E S was arbitrary and span(5a;) = H, it follows that AS = SA. 

Proposition 1: IfU is any unitary system for which W{U) ^ %, then for 
any xp e W{U) 

yV(U) = {Uip : U is a unitary operator in C^{U)} 

and the correspondence U -^ Uip is one-to-one. 

A Riesz basis for a Hilbert space H is the image under a bounded 
invertible operator of an orthonormal basis. Proposition 1 generalizes to 
generators of Riesz bases. A Riesz vector for a unitary system U is defined 
to be a vector tp for which Utp := {Utp : U E U} is a Riesz basis for the 
closed linear span of Uip, and it is called complete if span Uijj = H . Let 
TZWiU) denote the set of all complete Riesz vectors for U. 

Proposition 2: LetU he a unitary system on a Hilbert space H. If ip is a 
complete Riesz vector for U , then 

TZyV{U) ~ {A-ip : A is an operator in C^{U) that is invertible in B{H)}. 

2.4.3. Operator- Theoretic Interpolation 

Now suppose U is, a unitary system, such as IAd.t, and suppose 
{^iiip2T ■ ■ Tipm} C yV(Z^). (In the case of Ud,t, this means that 
("01: V'2, ■ • ■ 7 ipn) is au u-tuplc of wavelets. 

Let {Ai, A2, ■ ■ ■ , An) be an n-tuple of operators in the commutant U' of 
U, and let rj be the vector 



Then 



?7 := Aiipi + A21P2 H 1- Antpn ■ 



ry = AiVi + A^V^^'^i + ... AnV^y, 



(Ai+A2T/;;+... + A„F;;)Vi . (11) 
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We say that 77 is obtained by operator interpolation from 
{^pi,ip2, ■ ■ ■ ,1pm}- Since C^^{U) is a left W - module, it follows that the 
operator 

A:=A,+A,V^-^+...A,,v:};- (12) 

is an element of C.^^{U). Moreover, if B is another element of C^^iU) such 
that rj = B^i, then A- B e C^^{U) and {A - B)ipi = 0. So since Vi 
separates C^-^{U) it follows that A ^ B. Thus A is the unique element of 
C^-^{U) that takes ipi to rj. Let 5^i,...,^„ be the family of all finite sums of 
the form 

n 
i=0 

The is the left module of W generated by {/, V^^\ . . . , V^^}. It is the W- 
linear span of {/, V^^ , . . . , V^^}. Let 

■A^Vi,---,'/',! •= {S^u...,ii>Mi (13) 

So 

Al^i,...,^„ = <^A,V'^ : A^cU'\ . 

We call this the interpolation space for hi generated by {ipi, . . . , tpn)- From 
the above discussion, it follows that for every vector 77 e A^i/.i,i/.2,...,i/'„ there 
exists a unique operator A G C^^ (ZY) such that 77 = j47/;i , and moreover this 
A is an element of iS^i,...,^„. 

2.4.4. Normalizing the Commutant 

In certain essential cases (and we are not sure how general this type of 
case is) one can prove that an interpolation unitary V^ normalizes the 
commutant W of the system in the sense that V^U'Vjj = W . (Here, it is 
easily seen that [V^]* = V^ .) Wrfte V := V^. If V normahzes U\ then the 
algebra, before norm closure, generated by W and V is the set of all finite 
sums (trig polynomials) of the form ^ A„T^", with coefficients An S W , 
n ^ "L. The closure in the strong operator topology is a von Neumann 
algebra. Now suppose further that every power of V is contained in C^{U). 
This occurs only in special cases, yet it occurs frequently enough to yield 
some general methods. Then since C^iU) is a SOT-closed linear subspace 
which is closed under left multiplication by hi' , this von Neumann algebra 
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is contained in C^{IA)^ so its unitary group parameterizes a norm-path- 
connected subset of WiU) that contains i/; and -q via the correspondence 

In the special case of wavelets^ this is the basis for the work that Dai 
and I did in [11, Chapter 5] on operator-theoretic interpolation of wavelets. 
In fact, we specialized there and reserved the term operator-theoretic in- 
terpolation to refer explicitely to the case when the interpolation unitaries 
normalize the commutant. In some subsequent work, we loosened this re- 
striction yielding our more general definition given in this article, because 
there are cases of interest in which we weren't able to prove normalization. 
However, it turns out that if tp and 77 are s-elementary wavelets (see section 
2.5.4), then indeed V2 normalizes {D^T}' ■ (See Proposition 14.) Moreover, 
V2 has a very special form: after conjugating with the Fourier transform, 
it is a composition operator with symbol a natural and very computable 
measure-preserving transformation of M. In fact, it is precisely this special 
form for VV that allows us to make the computation that it normalizes 
{D, T}'. On the other hand, we know of no pair (-0, rf) of wavelets for which 
V2 fails to normalize {£',T}'. The difficulty is simply that in general it is 
very hard to do the computations. 

Problem: If {V', rj} is a pair of dyadic orthonormal wavelets, does the 
interpolation unitary V2 normalize {D, T}'? As mentioned above, the an- 
swer is yes if ip and rj are s-elementary wavelets. 



2.4.5. An Elementary Interpolation Result 

The following result is the most elementary case of operator-theoretic in- 
terpolation. 

Proposition 3: LetlA he a unitary system on a Hilbert space H. Ifipi and 
'02 are in WiU), then 

V'l + AV'2 e nw{u) 

for all complex scalars A with |A| 7^ 1. More generally, if tpi o,nd 1/12 are in 
TZWilA) then there are positive constants b > a > such that -01 -|- X'ip2 € 
TlWiU) for all X eC with either |A| < a or with |A| > b. 

Proof: If ipi,tp2 G W(Z^), let V be the unique unitary in C^^iU) given by 
Proposition 1 such that V^-)2 ~ ipi- Then 

■01 +A-02 = (y + A/)V'2- 
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Since V is unitary, (V + XI) is an invertiblc element of C^^iU) if |A| ^ 1, 
so the first conclusion follows from Proposition 2. Now assume i/'i,i/'2 G 
TZW{U). Let A be the unique invertible element oiC^p^i^) such that A'ip2 = 
Tpi, and write i/'i + -^V'2 = (^ + XI)'ip2- Since A is bounded and invertible 
there are 6 > a > such that 

aiA) C {z e C : a < |z| < 5} 

where (j{A) denotes the spectrum of A , and the same argument applies. D 

2.4.6. Interpolation Pairs of Wandering Vectors 

In some cases where a pair ip^rj oi vectors in yV{U) are given it turns out 
that the unitary V in C^{U) with Vip — rj happens to be a symmetry (i.e. 
V'^ = I). Such pairs are called interpolation pairs of wandering vectors, and 
in the case where U is a wavelet system, they are called interpolation pairs 
of wavelets. Interpolation pairs are more prevalent in the theory, and in 
particular the wavelet theory, than one might expect. In this case (and in 
more complex generalizations of this) certain linear combinations of com- 
plete wandering vectors are themselves complete wandering vectors - not 
simply complete Riesz vectors. 

Proposition 4: LetU he a unitary system, let tp,ri E W{U), and let V be 
the unique operator in C^ilA) with Vip — rj. Suppose 

V^ =/. 

Then 

cos a • ip + i sin a • 77 G WiJJ) 

for allO<a< 27r. 

The above result can be thought of as the prototype of our operator- 
theoretic interpolation results. It is the second most elementary case. More 
generally, the scalar a in Proposition 4 can be replaced with an appropriate 
self-adjoint operator in the commutant of U. In the wavelet case, after 
conjugating with the Fourier transform, which is a unitary operator, this 
means that a can be replaced with a wide class of nonnegative dilation- 
periodic (see definition below) bounded measurable functions on R. 
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2.4.7. A Test For Interpolation Pairs 

The following converse to Proposition 4 is typical of the type of computa- 
tions encountered in some wandering vector proofs. 

Proposition 5: Let hi be a unitary system, let iPtV ^ W(ZY), and let V be 
the unique unitary in C^{U) with Vip — rj. Suppose for some < a < ^ 
the vector 

p := cos a ■ ip + i sin a • rj 

is contained in W{U). Then 

Proof: Since Uip is a basis it will be enough to show that VUiip = V*Uiil} 
for all Ui E U. So it will suffice to prove that for all Ui,U2 gU we have 

Using the fact that V locally commutes with U at ip we have 

{VUi^, U2^) = {UiV^P, U2ik) = (C/ir/, t/2^) and 

{V*Uii>,U2-^) = {Ui^,VU2i^) - {Uii;,U2V^b) ^ {Ui^,U2v)- 

So we must show that {UxT], C/2V') — {Uii^, C^2'7) for all C/i, C/2 G U. 
Write p :— pa- By hypothesis xp^rj and p are unit vectors. So compute 

1 = {p,p) = cos^a • (V'jV') + i sin a cos a • (77, tp) 

i sin a cos a • (^, 77) + sin a ■ (rj, rj) 

= 1 + i sinacosa • ((?7, 7/;) - {ip^rf)). 

Thus, since sinacosa ^ 0, we must have (r/,')/') = (V':'?)- Also, for 
Ui,U2 gU with Ui 7^ U2 we have 

= (J7i/9, t/2/5) = cos^a • (J/ifA, f72'0) + isinacosa • (f7i?7, C/2V') 
- I sinacosa • (C/i-f/;, ^7277) + sin a ■ (Uirj , U2ri) 

~ isinacosa ■ {{Uirj,U2'(p) - (f/i-fA, f72??)), 
which implies {Uirj, C/2V') — {Ui'ip, C^2'7) as required. D 



February 2, 2008 5:34 WSPC/Lccturc Notes Scries: 9in x 6in larson-chaptcr 



Wavelet Sets and Operator- Theoretic Interpolation 17 

The above result gives an experimental method of checking whether 
V^ = I for a given pair ip,ri £ W{U). One just checks whether 

is an element of W(ZY), which is much simpler than attempting to work 
with the infinite matrix of V with respect to the basis Uip (or some other 
basis for H). 

2.4.8. Connectedness 

If we consider again the example of the left regular representation tt of a 
group G on H := 1^{G), then the local commutant oiU := 7r(G) at a vector 
ip € W{tt{G)) is just the commutant of ttIG). So since the unitary group 
of the von Neumann algebra (7r(G))' is norm-arcwise-connected, it follows 
that yV(7r(G')) is norm-arcwise-connected. 

Problem A in [11] asked whether W{D,L) is norm-arcwise-connected. It 
turned out that this conjecture was also formulated independently by Guido 
Weiss ([38], [37]) from a harmonic analysis point of view (our point of view 
was purely functional analysis), and this problem (and related problems) 
was the primary stimulation for the creation for the creation of the WUTAM 
CONSORTIUM - a team of 14 researchers based at Washington University 
and Texas A&M University. (See [52].) 

This connectedness conjecture was answered yes in [52] for the special 
case of the family of dyadic orthonormal MRA wavelets in L^ (M) , but still 
remains open for the family of arbitrary dyadic orthonormal wavelets in 
L^(R). 

In the wavelet case Ud,t , if "0 G W(D,T) then it turns out that 
C^{Ud,t) is in fact much larger than (Ud,tY = {D,Ty , underscoring 
the fact that Ud.t is NOT a group. In particular, {D,Ty is abelian while 
C^{'D,T) is nonabelian for every wavelet ip. (The proof of these facts are 
contained in [11].) 

2.5. Wavelet Sets 

2.5.1. The Fourier Transform 

We will use the following form of the Fourier-Plancherel transform JF on 
Ti = L^ (R) , because it is a form normalized so it is a unitary transformation. 
Although there is another such normalized form that is frequently used, and 
actually simpler, the present form is the one we used in our original first 
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paper [11] involving operator theory and wavelets, and so we will stick with 
it in these notes to avoid any confusion to a reader of both. 
If /,ge ii(M)ni2(K) then 

{Tf){s) := ^ / e— 7(t)rfi :- /(s), (14) 

and 

{T-'g){t) - ^ / e^''g{s)ds. (15) 

V ^TT Jv 



We have 

\/ l-K JR 

So TTaT-^g = e-^^^.g. For A e B{n) let i denote J-^J--!. Thus 

f„=Me-..., (16) 

where for h E L°° we use Mh to denote the multiplication operator / -^ hf. 
Since {M^-^o,s : a e R} generates the m.a.s.a. V{R) := {M,i : /i e L°°(R)} 
as a von Neumann algebra, we have 

Similarly, 

{TD''f){s) = -L / e-'"*(V2)"/(2"t)rfi 
v 27r jR 



V Stt Jr 
= (V2)-2(^/)(22-"^) = (D-^^Tf)(s). 

So 5" = D-" = D"\ Therefore, 

D^D^^^D*. (17) 

2.5.2. The Commutant of {D,T} 

We have T{D,T}'J=-'^ = {D,fy. It turns out that {D,fy has an easy 
characterization. 

Theorem 6: 

{D,fy = {Mh-. h e L°°(R) and /i(s) = /i(2s) a.c.}. 
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Proof: Since D = D* and D is unitary, it is clear that M^ G {D, T}' if 
and only if M^ commutes with D. So let g e L'^(R) be arbitrary. Then 
(a.e.) we have 

iMhDg){s) - h{s)iV2 g{2s)), and 
{DMHg){s) ^ Dihis)gis)) = Vh{2s)g{2s). 

Since these must be equal a.e. for arbitrary g, we must have h{s) — h{2s) 
a.e. D 

Now let E = [-2,-1) U [1, 2), and for n e Z let E„ = {2"x: x e E}. 
Observe that the sets £'„ are disjoint and have union E\{0}. So if g is 
any uniformly bounded function on E, then g extends uniquely (a.e.) to a 
function g E L°°(M.) satisfying 

.g(s) = g{2s), s € K, 

by setting 

5(2"s)=5(s), seE,neZ, 

and ^(0) — 0. We have H^Hoo — llslloo- Conversely, if h is any function 
satisfying h(s) — h{2s) a.e., then h is uniquely (a.e.) determined by its 
restriction to E. This 1-1 mapping g -^ Mg from L°°{E) onto {D^T}' is a 
^-isomorphism. 

We will refer to a function h satisfying h{s) = h(2s) a.e. as a 2-dilation 
periodic function. This gives a simple algorithm for computing a large class 
of wavelets from a given one, by simply modifying the phase: 

Given ip,, let tp = ^{ip), choose a real-valued function h G L°°(E) 
arbitrarily, let g — exp^ih), extend to a 2-dilation periodic (18) 

function g as above, and compute ipg = T^ (g^)- 

In the description above, the set E could clearly be replaced with 
[— 27r, — tt) U [tt, 27r), or with any other "dyadic" set [—2a, a) U [a, 2a) for 
some a > 0. 

2.5.3. Wavelets of Computationally Elementary Form 

We now give an account of s-elementary and M S'F-wavelets. The two most 
elementary dyadic orthonormal wavelets are the Haar wavelet and Shan- 
non's wavelet (also called the Littlewood-Paley wavelet). 
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The Haar wavelet is the function 



1, 0< t< i 
i^H{t) = { -1, i<i<l (19) 

0, otherwise. 

In this case it is very easy to sec that the dilates/translates 

{2tv-jj(2"-£): nJeZ} 

are orthonormal, and an elementary argument shows that their span is 
dense in L^{R). 

Shannon's wavelet is the L^(M)-function with Fourier transform ips — 
■^XEo where 

Eo = [-2tt, -tt) U [tt, 2tt). (20) 

The argument that Tps is a wavelet is in a way even more transparent than 
for the Haar wavelet. And it has the advantage of generalizing nicely. For 
a simple argument, start from the fact that the exponents 

{e*^": neZ} 

restricted to [0, 27r] and normalized by -7= is an orthonormal basis for 



L2[0, 2tt]. Write Eq = E_\JE+ where E^ = [-27r, -tt), E+ = [tt, 27r). Since 
{E- + 27r, Ej^} is a partition of [0, 27r) and since the exponentials e*^* are 
invariant under translation by 27r, it follows that 

: ne zl (21) 

Ea I 



/27r 
is an orthonormal basis for L'^{Eq). Since T — Mf,-is , this set can be written 

{f^^s- eez}. (22) 

Next, note that any "dyadic interval" of the form J = [b, 26), for some 6 > 
has the property that {2" J: n G Z}, is a partition of (0, 00). Similarly, any 
set of the form 

/C = [-2a, -a) U [6, 26) (23) 

for a, 6 > 0, has the property that 

{2"/C: neZ} 
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is a partition of R\{0}. It follows that the space L^(/C), considered as a 
subspace of L^ (R) , is a complete wandering subspace for the dilation unitary 
(£)/)(s) = V2 /(2s). For each n e Z, 

D"(L2(/C)) = i2(2-«/C). (24) 

So 0„£'"(i^(/C)) is a direct sum decomposition of L^{R). In particular 
Eq has this property. So 



L>" 



^ez^ = 



„2"ils 



2-"Bn 



^e z 



(25) 



is an orthonormal basis for L^(2 "i?o) for each n. It follows that 

is an orthonormal basis for L^(M). Hence {D'"'T^il;s- n,£ E Z} is an or- 
thonormal basis for L^(R), as required. 

The Haar wavelet can be generalized, and in fact Daubechie's well- 
known continuous compactly-supported wavelet is a generalization of the 
Haar wavelet. However, known generalization of the Haar wavelet are all 
more complicated and difficult to work with in hand-computations. 

For our work, in order to proceed with developing an operator algebraic 
theory that had a chance of directly impacting concrete function-theoretic 
wavelet theory we needed a large supply of examples of wavelets which 
were elementary enough to work with. First, we found another "Shannon- 
type" wavelet in the literature. This was the Journe wavelet, which we 
found described on p. 136 in Daubcchies book [14]. Its Fourier transform is 



Ej 



327r 



-47r U 



An 



U 



47r 
T 



■,7r U 



47r, ■ 



327r 



Then, thinking the old adage "where there's smoke there's fire!", we 
painstakingly worked out many more examples. So far, these are the basic 
building blocks in the concrete part of our theory. By this we mean the part 
of our theory that has had some type of direct impact on function-theoretic 
wavelet theory. 



2.5.4. Definition of Wavelet Set 

We define a wavelet set to be a measurable subset E of ] 



for which 



/2tt 



XE 



is the Fourier transform of a wavelet. The wavelet ipE '■— —^Xe is called 

V 27r 

s- elementary in [11]. 
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It turns out that this class of wavelets was also discovered and system- 
atically explored completely independently, and in about the same time 
period, by Guido Weiss (Washington University), his colleague and for- 
mer student E. Hernandez (U. Madrid), and his students X. Fang and X. 
Wang. In [17,37, 38] they are called MSF (minimally supported frequency) 
wavelets. In signal processing, the parameter s, which is the independent 
variable for -0, is the frequency variable, and the variable t, which is the 
independent variable for ip, is the time variable. No function with support 
a subset of a wavelet set E of strictly smaller measure can be the Fourier 
transform of a wavelet. 

Problem. Must the support of the Fourier transform of a wavelet contain 
a wavelet set? This question is open for dimension 1 . It makes sense for any 
finite dimension. 



2.5.5. The Spectral Set Condition 

^From the argument above describing why Shannon's wavelet is, indeed, a 
wavelet, it is clear that sufficient conditions for _E to be a wavelet set are 



(i) the normalized exponential ^=e , ^ G Z, when restricted 



1 ^Us 

to E should constitute an orthonormal basis for L'^{E) (in other 
words _E is a spectral set for the integer lattice Z), 



and 

(ii) The family {2"_E : n e Z} of dilates of E by integral powers of 
2 should constitute a measurable partition (i.e. a partition modulo 

nuU sets) of M. 

These conditions are also necessary. In fact if a set E satisfies (i) , then for 
it to be a wavelet set it is obvious that (ii) must be satisfied. To show that 
(i) must be satisfied by a wavelet set E, consider the vectors 



Since ipE is a wavelet these must be orthogonal, and so the sets 
{2"E: n E Z} must be disjoint modulo null sets. It follows that 
{^=e*^^|£;: £ e Z} is not only an orthonormal set of vectors in L'^{E), 



2ir 

2t 



it must also span L (E). 

It is known from the theory of spectral sets (as an elementary special 
case) that a measurable set E satisfies (i) if and only if it is a generator of a 
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measurable partition of M under translation by 27r (i.e. iff {E+2Trn : n G Z} 
is a measurable partition of M) . This result generalizes to spectral sets for 
the integral lattice in M". For this elementary special case a direct proof is 
not hard. 

2.5.6. Translation and Dilation Congruence 

We say that measurable sets E, F are translation congruent modulo 27r if 
there is a measurable bijection (p: E —^ F such that (/)(s) — s is an integral 
multiple of 27r for each s G _E; or cquivalcntly, if there is a measurable 
partition {_E„ : n G Z} of _E such that 

{£;„ + 2n7r : n G Z} (26) 

is a measurable partition of F. Analogously, define measurable sets G and 
H to be dilation congruent modulo 2 if there is a measurable bijection 
T : G ^ H such that for each s G G there is an integer n, depending on 
s, such that t{s) — 2"s; or cquivalcntly, if there is a measurable partition 
{Gnj'^oc of G such that 

{2"G}-^ (27) 

is a measurable partition of H. (Translation and dilation congruency mod- 
ulo other positive numbers of course make sense as well.) 
The following lemma is useful. 

Lemma 7: Let f G L^(R), and let E ~ supp(f). Then f has the property 
that 

{e'"'f: nGZ} 

is an orthonormal basis for L^ (E) if and only if 

(i) E is congruent to [0, 27r) modulo 27r, and 
(^^) l/(s)l = ^ a.e. onE. 

If £^ is a measurable set which is 27r-translation congruent to [0, 27r), 
then since 

'-= : £gZ 

/2lT [0,27r) 

is an orthonormal basis for L^[0, 27r] and the exponentials e*^'' are 27r- 
invariant, as in the case of Shannon's wavelet it follows that 

its 

£g 



E 
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is an orthonormal basis for L'^{E). Also, if E is 27r-translation congruent to 
[0, 27r), then since 

{[0,27r) + 27rn: n G Z} 

is a measurable partition of M, so is 

{E + 2Tm: ne Z}. 

These arguments can be reversed. 

We say that a measurable subset G C M is a 2- dilation generator of a 
partition of M if the sets 

2"G:={2"s: seG}, neZ (28) 

are disjoint and ]R\ U„ 2"G is a null set. Also, we say that ii^ C M is a 
2'K-translation generator of a partition of M if the sets 

E + 2m: ■.= {s + 2n'K: s<eE}, n E Z, (29) 

are disjoint and M\ U„ (E + 2n7r) is a null set. 

Lemma 8: A measurable set E <Z M is a 2'K-translation generator of a 
partition of M if and only if modulo a null set, E is translation congruent 
to [0, 27r) modulo 27r. Also, a measurable set G C M is a 2-dilation generator 
of a partition o/R if and only if modulo a null set, G is a dilation congruent 
modulo 2 to the set [—2tt, — tt) U [tt, 27r). 

2.5.7. A Criterion 

The following is a useful criterion for wavelet sets. It was published inde- 
pendently by Dai-Larson in [11] and by Fang-Wang in [17] at about the 
same time in December, 1994. In fact, it is amusing that the two papers 
had been submitted within two days of each other; only much later did we 
even learn of each others work and of this incredible timing. 

Proposition 9: Let E <Z R be a measurable set. Then E is a wavelet set 
if and only if E is both a 2-dilation generator of a partition (modulo null 
sets) of R and a 27: -translation generator of a partition (modulo null sets) 
of K. Equivalently, E is a wavelet set if and only if E is both translation 
congruent to [0, 27r) modulo 27r and dilation congruent to [— 27r, — 7r)U[7r, 27r) 
modulo 2. 



February 2, 2008 5:34 WSPC/Lecture Notes Series: 9in x 6in larson-chapter 



Wavelet Sets and Operator- Theoretic Interpolation 25 

Note that a set is 27r-translation congruent to [0, 2tt) iff it is 2tt- 
translation congruent to [— 27r,7r) U [tt, 27r). So the last sentence of Propo- 
sition]^ can be stated: A measurable set E is a. wavelet set if and only if 
it is both 27r-translation and 2-dilation congruent to the Littlewood-Paley 

set [-27r, -7r)U [tt, 27r). 

2.6. Phases 

If _E is a wavelet set, and if f{s) is any function with support E which 
has constant modulus — i= on E, then T~^{f) is a wavelet. Indeed, by 

Lemma [T] {T^ f : G Z} is an orthonormal basis for L^ [E) , and since the 
sets 2^E partition R, so L'^{E) is a complete wandering subspace for D, it 
follows that {D"T^f : n, ^ S Z} must be an orthonormal basis for L^(M), as 
required. In [17, 37, 38] the term MSF- wavelet includes this type of wavelet. 
So MSF-wavelets can have arbitrary phase and s-elementary wavelets have 
phase 0. Every phase is attainable in the sense of chapter 3 for an MSF or 
s-elementary wavelet. 

2.6.1. Some Examples of One- Dimensional Wavelet Sets 

It is usually easy to determine, using the dilation-translation criteria, in 
Proposition whether a given finite union of intervals is a wavelet set. In 
fact, to verify that a given "candidate" set £' is a wavelet set, it is clear 
from the above discussion and criteria that it suffices to do two things. 

(1) Show, by appropriate partitioning, that E is 2-dilation- 
congrucnt to a set of the form [—2a, —a) U [b, 2b) for some a,b > 0. 

and 

(2) Show, by appropriate partitioning, that E is 27r-translation- 
congruent to a set of the form [c, c -|- 2tt) for some real number 
c. 

On the other hand, wavelet sets suitable for testing hypotheses, can be 
quite difficult to construct. There are very few "recipes" for wavelet sets, 
as it were. Many families of such sets have been constructed for reasons 
including perspective, experimentation, testing hypotheses, etc., including 
perhaps the pure enjoyment of doing the computations - which are some- 
what "puzzle-like" in nature. In working with the theory it is nice (and 
in fact necessary) to have a large supply of wavelets on hand that permit 
relatively simple analysis. 
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For this reason we take the opportunity here to present for the reader 
a cohection of such sets, mainly taken from [11], leaving most of the "fun" 
in verifying that they are indeed wavelet sets to the reader. 

We refer the reader to [12] for a proof of the existence of wavelet sets 
in M^"' , and a proof that there are sufficiently many to generate the Borel 
structure of R^"^. These results are true for arbitrary expansive dilation 
factors. Some concrete examples in the plane were subsequently obtained by 
Soardi and Weiland, and others were obtained by Gu and Speegle. Two had 
also been obtained by Dai for inclusion in the revised concluding remarks 
section of our Memoir [11]. 

In these examples we will usually write intervals as half-open intervals 
[•, ) because it is easier to verify the translation and dilation congruency 
relations (1) and (2) above when wavelet sets are written thus, even though 
in actuality the relations need only hold modulo null sets. 

(i) As mentioned above, an example due to Journe of a wavelet which 
admits no multiresolution analysis is the s-elementary wavelet with wavelet 
set 



327r 



, -47r U 



An 



All 
y,7r|U 



47r, 



327r 



To see that this satisfies the criteria, label these intervals, in order, as 
Ji, J2, Js, J A and write J — UJi. Then 



Ji U AJ2 U 4 J3 U Ji 



327r IGtt 



7 



7 



U 



IGtt 327r 

— ,— 



This has the form [—2a, a) U [b, 26) so is a 2-dilation generator of a partition 
of M\{0}. Then also observe that 

{ Ji + 67r, J2 + 2tt, J3, J4 - Att} 

is a partition of [0, 27r). 

(ii) The Shannon (or Littlewood-Paley) set can be generalized. For any 
— TT < a < TT, the set 



S„ 



[-27r + 2a, -tt + a) U [tt + a, 27r + 2a) 



is a wavelet set. Indeed, it is clearly a 2-dilation generator of a partition of 
R\{0}, and to see that it satisfies the translation congruency criterion for 
— TT < a < (the case < a < tt is analogous) just observe that 



{[-2tt + 2a, 2tt) + Att, [-27r, -tt + a) + 2tt, [tt + a, 27r + 2a)} 
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is a partition of [0, 2tt). It is clear that ipE^ is then a continuous (in L^(I 
norm) path of s-elenientary wavelets. Note that 



lim Tpi 



1 



X127tA7t)- 



This is not the Fourier transform of a wavelet because the set [27r, 47r) is 
not a 2-dilation generator of a partition of R\{0}. So 

lim "tpEa 

is not an orthogonal wavelet. (It is what is known as a Hardy wavelet 
because it generates an orthonormal basis for _ff^(R) under dilation and 
translation.) This example demonstrates that W{D,T) is not closed in 

(iii) Journe's example above can be extended to a path. For — y < /? < 
J the set 



Jb 



327r 



-47r + 4/3 U 



-/?, 



47r\ 


'47r A 


47r + 4/3, 47r H 



is a wavelet set. The same argument in (i) establishes dilation congruency. 
For translation, the argument in (i) shows congruency to [4/3, 27r+4/3) which 
is in turn congruent to [0, 27r) as required. Observe that here, as opposed 
to in (ii) above, the limit of ^j^ as /3 approaches the boundary point j is 
a wavelet. Its wavelet set is a union of 3 disjoint intervals. 

(iv) Let A C [tt, ^) be an arbitrary measurable subset. Then there is 
a wavelet set W, such that W n [tt, ^) = A. For the construction, let 

B^ [27r, 37r)\2A, 



C 



-n,-^)\{A-2n) 



and D = 2A- An. 



Let 



W ^ 



Stt 



,2Tr]uAUB(JCUD. 



We have W D [tt, ^) = A. Observe that the sets [^, 27r), A, B, C, D, are 
disjoint. Also observe that the sets 

^,2ttYa,^B,2C,D, 
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are disjoint and have union [— 27r, — tt) U [tt, 27r). In addition, observe that 
the sets 

— ,2'k\,A,B-2t:,C + 27r, D + 27r, 

are disjoint and have union [0, 2tt). Hence VF is a wavelet set. 

(v) Wavelet sets for arbitrary (not necessarily integral) dilation factors 
other then 2 exist. For instance, if d > 2 is arbitrary, let 



A = 
B = 



2d7r 27r 



rf+1' d+1 
2tt 27r 



d? -V d+l 
2dTr 2d^TT 



d+l\P-l 



and \ct G = A\JB\JC. Then G is d- wavelet set. To see this, note that {A + 
2tt, B, C} is a partition of an interval of length 27r. So G is 27r-translation- 
congruent to [0, 2tt). Also, {A, B, d^^C} is a partition of the set [—da, — a)U 
[/?, d/3) for a = -;p^ , and /3 = ^jTrj : so from this form it follows that 
{d"G: n e Z} is a partition of M\{0}. Hence if V' := -^^H-^Xg), it 
follows that {d^^{d"'t — t): n,i £ Z} is orthonormal basis for L^(R), as 
required. 

2.7. Operator- Theoretic Interpolation of Wavelets: The 
Special Case of Wavelet Sets 

Let _B, F be a pair of wavelet sets. Then for (a.e.) x E E there is a unique 
y (z F such that x — y € 27rZ. This is the translation congruence property 
of wavelet sets. Also, for (a.e.) x E E there is a unique z £ F such that - is 
an integral power of 2. This is the dilation congruence property of wavelet 
sets. (See section 2.5.6.) 

There is a natural closed-form algorithm for the interpolation unitary 
VT'^ which maps the wavelet basis for 'tpE to the wavelet basis for 'tpp. 
Indeed, using both the translation and dilation congruence properties of 
{E,F}, one can explicitly compute a (unique) measure-preserving trans- 
formation (7 := a^ mapping R onto M. which has the property that VT^ is 
identical with the composition operator defined by: 

f^foa-' 
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for all / € L^(]R). With this formulation, compositions of the maps a be- 
tween different pairs of wavelet sets are not difhcult to compute, and thus 
products of the corresponding interpolation unitaries can be computed in 
terms of them. 



2.7.1. The Interpolation Map a 

Let E and F be arbitrary wavelet sets. Let a: E ^ F he the 1-1, onto map 
implementing the 27r-translation congruence. Since E and F both generated 
partitions of R\{0} under dilation by powers of 2, we may extend cr to a 
1-1 map of M. onto R by defining cr(0) = 0, and 

cr(s) = 2"cr(2-"s) for s G 2"£^, n € Z. (30) 

We adopt the notation a^ for this, and call it the interpolation map for the 
ordered pair (E, F). 

Lemma 10: In the above notation, a^ is a measure-preserving transfor- 
mation from R onto R. 

Proof: Let a :~ a^. Let i7 C R be a measurable set. Let ri„ = fi n T^E, 
n ^ Ij, and let En ~ 2^"il„ C E. Then {f2„} is a partition of fi, and 
wc have m{a{En)) = m{En) because the restriction of cr to £■ is measure- 
preserving. So 



m 



n n 

= J2 2"ni{a{En)) - ^ 2"m(i?„) 

n n 

= ^TO(2"S„)-^m(aO=™(0). ^ 

n n 

A function /: R ^ R is called 2-homogeneous if /(2s) — 2f{s) for all 
s e R. Equivalcntly, / is 2-homogeneous ifl: /(2"s) = 2"/(s), s € R, n e Z. 
Such a function is completely determined by its values on any subset of R 
which generates a partition of IR\{0} by 2-dilation. So a^ is the (unique) 2- 
homogeneous extension of the 27r-transition congruence E -^ F. The set of 
all 2-homogeneous measure-preserving transformations of R clearly forms 
a group under composition. Also, the composition of a 2-dilation-periodic 
function / with a 2-homogeneous function g is (in either order) 2-dilation 
periodic. We have /(ff(2s)) = fi2g{s)) = f{g{s)) and ff(/(2s)) = .g(/(s)). 
These facts will be useful. 
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2.7.2. An Algorithm For The Interpolation Unitary 

Now let 

UE--=U,-, (31) 

where if a is any measure-preserving transformation of R then Uo- denotes 
the composition operator defined by U„f = / o a^^, f E L^(K). Clearly 
(af )-i = ct| and (L/|)* = C/|. We have U0e = i^F since a§{E) = F. 
That is, 

Ue$e ^-ipEoaf^ -7k=^e ° erf = -7k=Xf = V-F. 
V^T'' v27r 

Proposition 11: Let E and F he arbitrary wavelet sets. Then U^ G 
Cz iD,T). Hence T^^U^T is the interpolation unitary for the ordered 
pair {-ipEjipp)- 

Proof: Write a ~ <t^ and Ua = U§. We have Uai's = V'-F since (j{E) — F. 
We must show 

U„Wf^$E = D'^f^U^^E, nJeZ. 
We have 

{U,D"f'$E)is) = iU„D"e-'''$E){s) 

= C/^2-^e-*'2-"'*^^(2-"s) 

This last term is nonzero iff 2 ^"s e f, in which case (T^^(2^"s) = (T^(2^"s) 
= 2~"s + 27rfc for some fc e Z since (t|; is a 27r-translation-congruence on 
F. It follows that e-'''^"'(2-"s) ^ e-^"-"s_ jjence we have 

{U„Wf^i^E){s) = 2-te-»''^"'">j.(2-"s) 
= (5"f '^j.)(s) 
= (5"f'C/,^i5)(s). 

We have shown U^ G C^ {D,T). Since U^'ipE = i^F, the uniqueness part 
of Proposition 1 shows that J^~^UE!F must be the interpolation unitary 
for {'4!e-,'4>f)- □ 
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2.8. The Interpolation Unitary Normalizes The 
Commutant 

Proposition 12: Let E and F he arbitrary wavelet sets. Then the inter- 
polation unitary for the ordered pair (ipE, V'-f) normalizes {D, T}' . 

Proof: By Proposition llll we may work with C/J)^ in tlic Fourier transform 
domain. By Theorem 6, the generic element of {D, T}' has the form Alh for 
some 2-dilation-periodic function h E L°°{M^). Write a = a^ and Ua = U^. 
Then 

U-HlhU, ^Ahoa-^. (32) 

So since the composition of a 2-dilation-periodic function with a 2- 
homogeneous function is 2-dilation-periodic, the proof is complete. D 

2.8.1. C^{D,T) is Nonabelian 

It can also be shown ([11, Theorem 5.2 (iii)]) that if E, F are wavelet sets 
with E ^ F then U^ is not contained in the double commutant {D, T}" . So 
since U^ and {D, T}' are both contained in the local commutant of ZYg f at 

tpE, this proves that Ct {D,T) is nonabelian. In fact (see [11, Proposition 
1.8]) this can be used to show that C^{D, T) is nonabelian for every wavelet 
V'. We suspected this, but we could not prove it until we discovered the 
"right" way of doing the needed computation using s-elementary wavelets. 
The above shows that a pair (S, F) of wavelets sets (or, rather, their cor- 
responding s-elementary wavelets) admits operator-theoretic interpolation 
if and only if GroupjC/^} is contained in the local commutant Cr {D,T), 
since the requirement that U^ normalizes {D,Ty is automatically satis- 
fied. It is easy to see that this is equivalent to the condition that for each 
n € Z, cr" is a 27r-congruence of E in the sense that (cr"'(s) — s)/2tt G Z for 
all s E E, which in turn implies that a"{E) is a wavelet set for all n. Here 
a = a^. This property hold trivially if a is involutive (i.e. a^ = identity). 

2.8.2. The Coefficient Criterion 

In cases where "torsion" is present, so (cr^)'^ is the identity map for some 
finite integer k, the von Neumann algebra generated by {Z?,T}' and U := 
U^ has the simple form 



,ra=0 



Mh,,U"': hn& L°°{R)withhn{2s) = hn{s), se 
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and so each member of this "interpolated" family of wavelets has the form 

k 



— J2hnis)Xa"{E) (33) 



27r „ 

n—0 

for 2-dilation periodic "coefficient" functions {hn{s)} which satisfy the nec- 
essary and sufhcient condition that the operator 

k 
^M;,„t/" (34) 

n=0 

is unitary. 

A standard computation shows that the map 6 sending ^q Mh^U" to 
the k X k function matrix (hij) given by 

hij = ^a(ij) ° 0-"'+^ (35) 

where a{i,j) — (J + 1) modulo fc, is a *-isomorphisin. This matricial algebra 
is the cross-product of {D, T}' by the *-automorphism ad{U^) correspond- 
ing to conjugation with U^. For instance, if fc = 3 then 9 maps 

Mh,+Mh,Ul:+Mh,{U^f 
to 

hi h2 ha \ 

/i3 o a^^ hi o a^^ /i2 o <t^^ . (36) 

y/l2 O (T^^ h^ o a^^ hi o <T^^ J 

This shows that ^g MhJJ"^ is a unitary operator iff the scalar matrix 
{hij){s) is unitary for almost all s S R. Unitarity of this matrix- valued 
function is called the Coefficient Criterion in [11], and the functions hi are 
called the interpolation coefficients. This leads to formulas for families of 
wavelets which are new to wavelet theory. 

2.9. Interpolation Pairs of Wavelet Sets 

For many interesting cases of note, the interpolation map a^ will in fact be 
an involution of M (i.e. aoa = id, where a :— a^, and where id denotes the 
identity map). So torsion will be present, as in the above section, and it will 
be present in an essentially simple form. The corresponding interpolation 
unitary will be a symmetry in this case (i.e. a selfadjoint unitary operator 
with square /). 

It is curious to note that verifying a simple operator equation U^ — I 
directly by matricial computation can be extremely difficult. It is much 
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more computationally feasible to verify an equation such as this by point- 
wise (a.e.) verifying explicitly the relation a o a — id for the interpolation 
map. In [11] we gave a number of examples of interpolation pairs of wavelet 
sets. We give below a collection of examples that has not been previously 
published: Every pair sets from the Journe family is an interpolation pair. 

2.10. Journe Family Interpolation Pairs 

Consider the parameterized path of generalized Journe wavelet sets given 
in [11, Example 4.5(iii)]. We have 



J 



P 



-^,-47r-4/3)U 



47r\ 



47r 

y,7r + /3|U 



47r + 4/3, 47r , 

7 



where the set of parameters j3 ranges — y < /3 < y . 
Proposition 13: Every pair (J^-^^Jfi^) is an interpolation pair. 

Proof: Let /3i,/32 G [-f , y) with /3i < /32. Write a = a^f^^ . We need to 
show that 

a2(x)=x (*) 

for all a; € M. Since a is 2-homogeneous, it suffices to verify (*) only for 
X £ J/3i . For X G Jp^ C\ Jp^ we have (j{x) = x, hence cr^(a;) = x. So we only 
need to check (*) for x e {J i3x\J p-i) ■ We have 

JpAJp2 = [-7r + /3i,-7r + /32)U[47r + 4/3i,47r + 4/32). 

It is useful to also write 

Jr32\Jr3i - h47r + 4/3i, -At: + Ap^) U [^ + /3i, tt + P^)- 

On [— 7r+/3i, — 7r+/32) we have (t(x) = a;+27r, which lies in [7r+/3i, 7r+/32). 
If we multiply this by 4, we obtain 4cr(x) e [At: + A(3i,Att + 4/32) C Jp^. 
And on [Att + 4/3i, 47r + 4/32) we clearly have a-{x) = x — 8tt, which lies in 

[-47r + 4/3i,-47r + 4/32). 

So for X e [— TT + /3i, — TT + (32) we have 

cr^(x) = a{a{x)) = -cr(4cr(a;)) = -[4cr(x)-87r] = cr(x)-27r = x+2n-2n = x. 

On [47r + 4/3i, 47r + 4/32) we have cr(x) = a; — Stt, which lies in [—Att + 
4/3i, -Att + 4/32). So icr(x) G [-tt + /3i, -tt + [32). Hence 

c^ (^cr(x)j ^ -<t{x) + 2tt 
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and thus 



-a{x) + 2tt 
4 



(j{x) + Stt — X ~ 8tt + Stt = X 



as required. 

We have shown that for all x E Jp^ we have cr'^{x) = x. This proves that 
(J^j, J^j) is an interpolation pair. D 

3. Unitary Systems and Frames 

In [33] we developed an operator-theoretic approach to discrete frame the- 
ory (i.e. frame sequences, as opposed to continuous frame transforms) on 
a separable Hilbert space. We then applied it to an investigation of frame 
vectors for unitary systems, frame wavelets and group representations. The 
starting-point idea, which is pretty simple-minded in fact, is to realize any 
frame sequence for a Hilbert space H as a compression of a Riesz basis for 
a larger Hilbert space. In other words, a frame is a sequence of vectors in a 
Hilbert space which dilates, (in the operator-theoretic or geometric sense, 
as opposed to the function-theoretic sense of multiplication of the indepen- 
dent variable of a function by a dilation constant), or extends, to a (Riesz) 
basis for a larger space. From this idea much can be developed, and some 
new perspective can be given to certain concepts that have been used in 
engineering circles for many years. See section l!^.2. 1 below. 

3.1. Basics on Frames 

Let H he a, separable complex Hilbert space. Let B{H) denote the algebra 
of all bounded linear operators on H . Let N denote the natural numbers, 
and Z the integers. We will use JJ to denote a generic countable (or finite) 
index set such as Z, N, Z^^) , N U N etc. 

A sequence {xj : j e N} of vectors in H is called a frame if there are 
constants A,B > {) such that 



A\\xr<y\{x,x,)\'<B\\x\\ 



for all X E H . The optimal constant (maximal for A and minimal for B) 
are called the /rame bounds. The frame {xj} is called a tight frame if 
A — B, and is called Parseval \i A = B = 1. (Originally, in [33] and a in 
number of subsequent papers, the term normalized tight frame was used 
for this. However, this term had also been applied by Benedetto and Ficus 
[5] for another concept: a tight frame of unit vectors; what we now call a 
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uniform tight frame, or spherical frame. So, after all parties involved, the 
name Parseval was adopted. It makes a lot of sense, because a Parseval 
frame is precisely a frame which satisfies Parseval's identity.) A sequence 
{xj} is defined to be a Riesz basis if it is a frame and is also a basis for H 
in the sense that for each x € H there is a unique sequence {aj} in C such 
that a; = ^ Oij^j with the convergence being in norm. We note that a Riesz 
basis is also defined to be basis which is obtained from an orthonormal basis 
by applying a bounded linear invertible operator. This is equivalent to the 
first definition. It should be noted that in Hilbert spaces the Riesz bases 
are precisely the bounded unconditional bases. We will say that frames 
{xj : j e J} and {j/j : j € JJ} on Hilbert spaces H, K, respectively, are 
unitarily equivalent if there is a unitary U : H ^ K such that Uxj = yj 
for all j € J. We will say that they are similar (or isomorphic) if there is a 
bounded linear invertible operator T: H ^ K such that Txj = yj for all 

Example 14: Let K = _L^(T) where T is the unit circle and measure is 
normalized Lebesgue measure, and let {e*"^ : n E Z} he the standard 
orthonormal basis for L^(T). If i<^ C T is any measurable subset then 
{6"'"!^: n e Z} is a Parseval frame for L'^{E). This can be viewed as 
obtained from the single vector Xe ^Y ^PPlyii^g s-U integral powers of the 
(unitary) multiplication operator Mgis . It turns out that these are all (for 
different E) unitarily inequivalent. This is an example of a Parseval frame 
which is generated by the action of a unitary group on a single vector. This 
can be compared with the definition of a frame wavelet. (As one might ex- 
pect, a single function ip in L^ (R) which generates a frame for L^ (R) under 
the action ofUo^T is called a frame- wavelet .j 



3.2. Dilation of Frames: The Discrete Version of 
Naimark's Theorem 

Now let {x„}„gj be a Parseval frame and let 0: H ^ K :— P{1) be 
the usual analysis operator (this was called the frame transform in [HL]) 
defined by 0{x) := ((x, a;„))„gj. This is obviously an isometry. Let P be the 
orthogonal projection from K onto 0(H). Denote the standard orthonormal 
basis for P{I) by {cj : j G J}. For any to e J, we have 

{6{x„i),Pen) = {P6{xm),en) = (6'(x,„),e„) 



February 2, 2008 5:34 WSPC/Lccturc Notes Scries: 9in x 6in larson-chaptcr 



36 Larson 

It follows easily that 0(x„) — Pcn, n E J. Identifying H with 0{H), this 
shows indeed that every Parseval frame can be realized by compressing an 
orthonornial basis, as claimed earlier. 

This can actually be viewed as a special case (probably the simplest 
possible special case) of an old theorem of Naimark concerning operator 
algebras and dilation of positive operator valued measures to projection 
valued measures. The connection between Naimark's theorem and the di- 
lation result for Parseval frames, and that the latter can be viewed as a 
special case of the former, was pointed out to me by Chandler Davis and 
Dick Kadison in a conference (COSY-1999: The Canadian Operator Alge- 
bra Symposium, Prince Edward Island, May 1999). 

3.3. Complements of Frames 

It is useful to note that P will equal I iff {x„} is a basis. Indeed, if P 7^ /, 
then choose 2; 7^ 0, z e (/ — P)K, and write z = ^ OinS-n for some sequence 
a„ G C. Then = Pz — ^q:„0(x„), and not all the scalars a„ are zero. 
Hence {x„} is not topologically linearly independent so cannot even be a 
Schauder basis. On the other hand if P = / then {x„} is obviously an 
orthonormal basis. 

Suppose {a;„}„gj is a Parseval frame for iJ, and let 9,P,K,en be as 
above. Let M = (I — P)K. Then j/„ :— (I — P)en is a Parseval frame on 
M which is complementary to {a;„} in the sense that the inner direct sum 
{Xn (B iJn' ri e J} is an orthonormal basis for the direct sum Hilbert space 
H © M. Moreover there is uniqueness: The extension of a tight frame to an 
orthonormal basis described in the above paragraph is unique up to unitary 
equivalence. That is if N is another Hilbert space and {zn} is a tight frame 
for N such that {a;„ ® Zn : n G J} is an orthonormal basis for H (B N, then 
there is a unitary transformation U mapping M onto N such that Uyn ~ z„ 
for all n. In particular, diniM = dim A^. 

If {xj} is a Parseval frame, we will call any Parseval frame {zj} such 
that {xj © Zj} is an orthonormal basis for the direct sum space, a strong 
complement to {xj}. So every Parseval frame has a strong complement 
which is unique up to unitary equivalence. More generally, if {yj} is a 
general frame we will call any frame {wj} such that {yj © Wj} is a Riesz 
basis for the direct sum space a complementary frame (or complement) to 

The notion of strong complement has a natural generalization. Let 
{x„}„gj and {ynjneJ be Parseval frames in Hilbert spaces H,K, respec- 
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tively, indexed by the same set J. Call these two frames strongly disjoint 
if the (inner) direct sum {xn ® Vn'- n e J} is a Parseval frame for the 
direct sum Hilbert space H (B K. It is not hard to see that this property 
of strong disjointness is equivalent to the property that the ranges of their 
analysis operators are orthogonal in P(J). More generally, we call a fc-tuple 
of Parseval frames ({zi„}„gj, . . . , {zknjnes) in Hilbert spaces Hi,. . . ,Hk, 
respectively, a strongly disjoint fc-tuple if {zi„ • • • ® Zkn '■ n G 1} is 
a Parseval frame for Hi © • • • © Hk , and we call it a complete strongly 
disjoint fc-tuple if {zi„ • • • © Zkn '■ n g JJ} is an orthonormal basis for 
Hi (B ■ ■ ■ ® Hk- li Ot: Hi -^ P{^) is the frame transform, 1 < i < fc, then 
strong disjointness of a fc-tuple is equivalent to mutual orthogonality of 
{ran 9i: 1 < i < fc}, and complete strong disjointness is equivalent to the 

k 

condition that ran 9i — P{I). 

There is a particularly simple intrinsic (i.e. non-geometric) characteri- 
zation of strong disjointness which is potentially useful in applications: Let 
{a;„}„gj and {yn}nei be Parseval frames for Hilbert spaces H and K, re- 
spectively. Then {x„} and {y-n} are strongly disjoint if and only if one of 
the equations 

V^(a;, Xn)yn = for all x e iJ (37) 

nea 

or ^{y^ yn)x„ = for all y e K 
holds. Moreover, if one holds the other holds also. 

3.4. Super-frames, Super-wavelets, and Multiplexing 

Suppose that {x„}„gj and {yn}nel are strongly disjoint Parseval frames 
for Hilbert spaces H and K, respectively. Then given any pair of vectors 
X & H, y ^ K, we have that 

X^^{x,Xn)Xn, V = ^{y,yn)yn- 

n n 

If we let a„ = {x, Xn) and &„ == {y, j/„), and then let c„ = a„ + bn, we have 

y^ a«2/n =0, ^ bnXn = 0, 

n n 

by (|37(l and therefore we have 

X = ^c„Xn, y = ^c„y„. (38) 

n N 
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This says that, by using one set of data {c„}, we can recover two vectors 
X and y (they may even he in different Hilbert spaces) by applying the 
respective inverse transforms (synthesis operators) corresponding to the 
two frame {xn} and {yn}- The above argument obviously extends to the 
fc-tuple case: If {fin- n e J}, i = 1, . . . ,k, is a strongly disjoint fc-tuple 
of Parseval frames for Hilbert spaces Hi, . . . ,Hk, and if {xi, . . . ,Xk) is an 
arbitrary fc-tuple of vectors with Xi ^ Hi, 1 < i < k, then 1)38(1 generalizes 
to 

•^2 — / ^ x-^i^ Jin) Jin 

neJI 

for each 1 < i < fc. So if we define a single "master" sequence of complex 
numbers {c„ : n e JJ} by 

fe 

Cn ^ / ^ X^ii Jin) ; 
i=l 

then the strong disjointness implies that for each individual i we have 



neJI 



C-nJii 



This simple observation might be useful in applications to data compression. 

In [33] we called such an n-tuple of strongly disjoint (or simply just 
disjoint) frames a super-frame, because it (or rather its inner direct sum) is 
a frame for the superspace which s the direct sum of the individual Hilbert 
spaces for the frames. In connection with wavelet systems this observation 
lead us to the notion of superwavelet, which is a particular type of vector- 
valued wavelet. In operator-theoretic terms this is just a restatement of the 
fact outlined above that a strongly disjoint fc-tuple of Parseval frames have 
frame-transforms which are isometries into the same space P(JJ) which have 
mutually orthogonal ranges. 

The notion of superframes and superwavelets, and many of their prop- 
erties, were also discovered and investigated by Radu Balan [3] in his Ph.D. 
thesis, in work that was completely independent from ours. 

3.5. Frame Vectors For Unitary Systems 

Let U he a. unitary system on a Hilbert space H. Suppose W{U) is 
nonempty, and fix V' G W{U). Recall from Section 1 that if t] is an arbi- 
trary vector in H, then rj G W{U) if and only if there is a unitary V (which 
is unique if it exists) in the local commutant C^iJJ) such that Vip = rj. 
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The following proposition shows that this idea generalizes to the theory of 
frames. Analogously to the notion of a wandering vector and a complete 
wandering vector, a vector x Cz H is called a Parseval frame vector (resp. 
frame vector with bounds a and b) for a unitary system hi if Ux forms a 
tight frame (resp. frame with bounds a and b) for spanilAx). It is called a 
complete Parseval frame vector (resp. complete frame vector with bounds 
a and h) when lAx is a Parseval frame (resp. frame with bounds a and b) 
for H. 

Proposition 15: Suppose that tjj is a com,plete wandering vector for a 
unitary system hi. Then 

(i) a vector rj is a Parseval frame vector for U if and only if there is a 

(unique) partial isometry A € C^(U) such that Aip ~ r/. 
(ii) a vector rj is a complete Parseval frame vector for hi if and only if 
there is a (unique) co-isometry A £ C^ihi) such that Aijj = rj. 

The above result does not tell the whole story. The reason is that many 
unitary systems do not have wandering vectors but do have frame vectors. 
For instance, this is the case in Example II 41 where the unitary system is the 
group of multiplication operators hi = {M^ins : n E Z,} acting on L'^{E). In 
the case of a unitary system such as the wavelet system hiD,T there exist 
both complete wandering vectors and nontrivial Parseval frame vectors, so 
the theory seems richer (however less tractable) and Proposition II 51 is very 
relevant. 

Much of Example El generalizes to the case of an arbitrary countable 
unitary group. There is a corresponding (geometric) dilation result. 

Proposition 16: Suppose that U is a unitary group such that WilA) is 
non-empty. Then every complete Parseval frame vector must be a complete 
wandering vector. 

Theorem 17: Suppose that U \s a. unitary group on H and 77 is a complete 
Parseval frame vector for U. Then there exists a Hilbert space K Z) H and 
a unitary group Q on K such that Q has complete wandering vectors, H is 
an invariant subspace of Q such that Q\h = hi, and the map g -^ g\H is a 
group isomorphism from Q onto hi. 

The following is not hard, but it is very useful. 

Proposition 18: Suppose that hi is a unitary group which has a complete 
Parseval frame vector. Then the von Neumann algebra w*{hl) generated by 
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U is finite. 

3.6. An Operator Model 

The following is a corollary of Theorem 1 171 It shows that Example 1141 can 
be viewed as a model for certain operators. 

Corollary 19: Let T G B{H) he a unitary operator and let rj ^ H be a 
vector such that {T^iy. n ^ "L} is a Parseval frame for H . Then there is a 
unique (modulo a null set) measurable set E <zT such that {T^rj: n G Z} 
and {e™'*!^; : n G Z} are unitarily equivalent frames. 

3.7. Group Representations 

These concepts generalize. For a unitary system hi on a Hilbert space iJ, 
a closed subspace M of iJ is called a complete wandering subspace for U 
if span {UM: U € U} is dense in H, and UM _L VM with U ^V. Let 
{si: i G /} be an orthonormal basis for M. Then AI is a complete wan- 
dering subspace for U if and only if {Uei : U EU,i E 1} is an orthonormal 
basis for H. We call {ci} a complete multi-wandering vector. Analogously, 
an n-tuple (ryi, . . . , ry„) of non-zero vectors (here n can be cxd) is called com- 
plete Parseval multi-frame vector for U if {Urji : U E U, i = 1, . . . ,n} forms 
a complete Parseval frame for H. Let G be a group and let A be the left 
regular representation of G on P'{G). Then {Ag x /„ : 5 G G} has a complete 
multi-wandering vector (/i,...,/„), where /i = (xg, 0, . . . ,0), . . . , /„ = 
(0, 0, ... , Xe). Let P be any projection in the commutant of (A ® In){G). 
Then (P/i, . . . , Pfn) is a complete Parseval multi-frame vector for the sub- 
representation (A (8) In)\p- It turns out that every representation with a 
complete Parseval multi-frame vector arises in this way. Item (i) of the fol- 
lowing theorem is elementary and was mentioned earlier; it is included for 
completeness. 

Theorem 20: Let G be a countable group and let tt be a representation 
of G on a Hilbert space H. Let A denote the left regular representation of 
G on P{G). Then 

(i) if 7r(G) has a complete wandering vector then n is unitarily equiv- 
alent to A, 

(ii) if 7r(G) has a complete Parseval frame vector then n is unitarily 
equivalent to a subrepresentation of A, 

(iii) if 7r(G) has a complete Parseval multi-frame vector 
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{V'l) 4'2, ■ ■ ■ ,4'n}j for some 1 < n < oo, then tt is unitarily equiva- 
lent to a subrepresentation of A (g) /„. 

4. Decompositions of Operators and Operator- Valued 
Frames 

The material we present here is contained in two recent papers. The first 
[f 5] was authored by a [VfGRE/REU] team consisting of K. Dykema, D. 
Freeman, K. Kornelson, D. Larson, M. Ordower, and E. Weber, with the 
title Ellipsoidal Tight Frames. This article started as an undergraduate re- 
search project at Texas A&M in the summer of 2002, in which Dan Freeman 
was the student and the other five were faculty mentors. Freeman is now 
a graduate student at Texas A&M. The project began as a solution of a 
finite dimensional frame research problem, but developed into a rather tech- 
nically deep theory concerning a class of frames on an infinite dimensional 
Hilbert space. The second paper [44], entitled Rank-one decomposition of 
operators and construction of frames, is a joint article by K. Kornelson and 
D. Larson. 

4.1. Ellipsoidal Frames 

We will use the term spherical frame (or uniform frame) for a frame se- 
quence which is uniform in the sense that all its vectors have the same 
norm. Spherical frames which are tight have been the focus of several ar- 
ticles by different researchers. Since frame theory is essentially geometric 
in nature, from a purely mathematical point of view it is natural to ask: 
Which other surfaces in a finite or infinite dimensional Hilbert space con- 
tain tight frames? (These problems can make darn good REU projects, in 
particular.) In the first article we considered ellipsoidal surfaces. 

By an ellipsoidal surface we mean the image of the unit sphere 5*1 in the 
underlying Hilbert space H under a bounded invertible operator A in B(H), 
the set of all bounded linear operators on H. Let Ea denote the ellipsoidal 
surface Ea ■— ASi. A frame contained in Ea is called an ellipsoidal frame, 
and if it is tight it is called an ellipsoidal tight frame (ETF) for that surface. 
We say that a frame bound K is attainable for Ea if there is an ETF for 
Ea with frame bound K. 

Given an ellipsoidal surface E := Ea, we can assume E = Et where 
T is a positive invertible operator. Indeed, given an invertible operator A, 
let A* = U\A*\ be the polar decomposition, where |^*| = (AA*)''-^^. Then 
A= \A*\U*. By taking T = |A*|, we see that TSi = ASi. Moreover, it is 
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easily seen that the positive operator T for which E = Ej- is unique. 

The starting point for the work in the first paper was the following 
Proposition. For his REU project Freeman found an elementary calculus 
proof of this for the real case. Others have also independently found this 
result, including V. Paulsen, and P. Casazza and M. Leon. 

Proposition 21: Let Ea be an ellipsoidal surface on a finite dimensional 
real or complex Hilbert space H of dimension n. Then for any integer k> n, 
Ea contains a tight frame of length k, and every ETF on Ea of length k 
has frame hound K = k [irace(r^^)l 

We use the following standard definition: For an operator B E H, the 
essential norm of B is: 

||i?||ess := inf{||_B — K\\ : K is a compact operator mB{H)} 

Our main frame theorem from the first paper is: 

Theorem 22: Let Ea be an ellipsoidal surface in an infinite dimensional 
real or complex Hilbert space. Then for any constant K > ||r~^||~^, Et 
contains a tight frame with frame bound K . 

So, for fixed A, in finite dimensions the set of attainable ETF frame 
bounds is finite, whereas in infinite dimensions it is a continuum. 

Problem. If the essential norm of A is replaced with the norm of A in the 
above theorem, or if the inequality is replaced with equality, then except for 
some special cases, and trivial cases, no theorems of any degree of generality 
are known concerning the set of attainable frame bounds for ETF's on Ea- 
It would be interesting to have a general analysis of the case where A — I 
is compact. In this case, one would want to know necessary and sufficient 
conditions for existence of a tight frame on Ea with frame bound 1. In the 
special case A = I then, of course, any orthonormal basis will do, and these 
are the only tight frames on Ea in this case. What happens in general when 
II A|| ess = 1 and A is a small perturbation of /? 

We use elementary tensor notation for a rank-one operator on H. Given 
u,v,x Cz H, the operator m (g) w is defined by {u (g) v)x — {x, v)u for x Cz H. 
The operator u u is a projection if and only if \\u\\ = 1. 

Let {xj}j be a frame for H. The standard frame operator is defined by: 
Sw — '^j{w,Xj)xj — ^- (xj (^Xj)w . Thus S — J2j^j '^ ^ji where this 
series of positive rank-1 operators converges in the strong operator topology 
(i.e. the topology of pointwise convergence). In the special case where each 
Xj is a unit vector, S is the sum of the rank-1 projections Pj ~ Xj Xj. 
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For A a positive operator, we say that A has a projection decomposition 
if A can be expressed as the sum of a finite or infinite sequence of (not 
necessarily mutually orthogonal) self-adjoint projections, with convergence 
in the strong operator topology. 

If Xj is a frame of unit vectors, then S — "^jXj (^ Xj is a projection 
decomposition of the frame operator. This argument is trivially reversible, 
so a positive invertible operator S is the frame operator for a frame of unit 
vectors if and only if it admits a projection decomposition S = ^ Pj. If the 
projections in the decomposition are not of rank one, each projection can 
be further decomposed (orthogonally) into rank-1 projections, as needed, 
expressing S = ^^ x„ (g) x„, and then the sequence {x„} is a frame of unit 
vectors with frame operator S. 

In order to prove Theorem 22, we first proved Theorem 23 (below), 
using purely operator-theoretic techniques. 

Theorem 23: Let A be a positive operator in B{H) for H a real or complex 
Hilbert space with infinite dimension, and suppose ||j4||es5 > 1. Then A has 
a projection decomposition. 

Suppose, then, that {x„} is a frame of unit vectors with frame operator 
S. If we let j/j = S~^Xj, then {yj}j is a Parseval frame. So {yj}j is an 
ellipsoidal tight frame for the ellipsoidal surface E _i ~ S^'^Si. This ar- 
gument is reversible: Given a positive invertible operator T, let S = T^^. 
Scale T if necessary so that ||S'||ess > 1- Let S = ^ Xj®Xj be a projection 
decomposition of S. Then {Txj} is an ETF for the ellipsoidal surface TSi. 
Consideration of frame bounds and scale factors then yields Theorem 22. 

Most of our second paper concerned weighted projection decomposi- 
tions of positive operators, and resultant theorems concerning frames. If T 
is a positive operator, and if {c„} is a sequence of positive scalars, then a 
weighted projection decomposition of T with weights {c„} is a decompo- 
sition T = ^ Pj where the Pj are projections, and the series converges 
strongly. We have since adopted the term targeted to refer to such a de- 
composition, and generalizations thereof. By a targeted decomposition of T 
we mean any strongly convergent decomposition T — ^^ r„ where the r„ 
is a sequence of simpler positive operators with special prescribed proper- 
ties. So a weighted decomposition is a targeted decomposition for which the 
scalar weights are the prescribed properties. And, of course, a projection 
decomposition is a special case of targeted decomposition. 

After a sequence of Lemmas, building up from finite dimensions and 
employing spectral theory for operators, we arrived at the following theo- 
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rem. We will not discuss the details here because of limited space. It is the 
weighted analogue of theorem 23. 

Theorem 24: Let i? be a positive operator in B{H) for H with ||_B||ess > 
1. Let {ci\°°^i be any sequence of numbers with < q < 1 such that 
^j Ci = oo. Then there exists a sequence of rank-one projections {Pi\^i 
such that B =- ^Y^=\ '^iPi 



4.2. A Problem in Operator Theory 

We will discuss a problem in operator theory that was motivated by a prob- 
lem in the theory of Modulation Spaces. We tried to obtain an actual "re- 
formulation" of the modulation space problem in terms of operator theory, 
and it is well possible that such a reformulation can be found. At the least 
we (Chris Heil and myself) found the following operator theory problem, 
whose solution could conceivably impact mathematics beyond operator the- 
ory. I find it rather fascinating. I need to note that we subsequently showed 
(in an unpublished jointly-written expository article) that the actual modu- 
lation space connection requires a modified and more sophisticated version 
of the problem we present below. I still feel, that the problem I will present 
here has some independent interest, and may serve as a "first step" in de- 
veloping a theory that might have some usefulness. Thus, I hope that the 
reader will find it interesting. 

Let H be an infinite dimensional separable Hilbert space. As usual, 
denote the Hilbert space norm on H by || • ||. If x and y are vectors in H, 
then xi^y will denote the operator of rank one defined by {xiSiy)z = (z, y)x. 
The operator norm of a; (g) y is then just the product of ||a;|| and \\y\\. 

Fix an orthonormal basis {e„}„ for H. For each vector v in H, define 

\M\^Y.\^v,e„)\ 

n 

This may be -|-cx3. 

Let L be the set of all vectors w in H for which || |v|| | is finite. Then L is 
a dense linear subspace of H, and is a Banach space in the "triple norm". 
It is of course isomorphic to £^ 

Let T be any positive trace-class operator in B{H). 
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The usual eigenvector decomposition for T expresses T as a series con- 
verging in the strong operator topology of operators hn ® hm where {/in} 
is an orthogonal sequence of eigenvectors of T . That is, 



T = ^K®hn 



In this representation the eigenvalue corresponding to the eigenvector ft,„ is 
the square of the norm: ||/i„|p. The trace of T is then 



Eii'^" 



and since T is positive this is also the trace-class norm of T . 

Let us say that T is of Type A with respect to the orthonormal basis 
{e„} if, for the eigenvectors {/i„} as above, we have that X^nlll^nllP i^ 
finite. [Note that this is just the (somewhat unusual) formula displayed 
above for the trace of T with the triple norm used in place of the usual 
Hilbert space norm of the vectors {ft.„}.] 

And let us say that T is of Type B with respect to the orthonormal 
basis {e„} if there is some sequence of vectors {w„} in H with ^^ IlknllP 
finite such that 

T^^Vn(S)hn 
n 

where the convergence of this series is in the strong operator topology. 

Problem: If T if of Type B with respect to an orthonormal basis {e„}, 
then must it be of Type A with respect to {e„}? 

Note: If the answer to this problem is negative (as I suspect it is), then 
the following subproblcm would be an interesting one. 

Subproblem: Let {e„} be an orthonormal basis for H. Find a char- 
acterization of all positive trace class operators T that are of Type B with 
respect to {e„}. In particular, is every positive trace class operator T of 
Type B with respect to {e„}? My feeling is no. (See the next example.) 

Example 25: Let x be any vector in H that is not in L, and \etT — x®x. 
Then T is trace class, in fact has rank one, but clearly T is clearly not of 
Type A. Can such a T be of type B? (I don't think it is necessarily of Type 
B for all such T, however.) 
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